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Abstract. We give an asymptotic expansion {the higher Stirling 
formula) and an infinite product representation (the Weierstrass 
product representation) of tlie Vigneras multiple gamma functions 
by considering the classical limit of the multiple g-gamma func- 
tions. 



1. Introduction 

The multiple gamma function was introduced by Barnes. It is de- 
fined to be an infinite product regularized by the multiple Hurwitz zeta 
functions [|], [|], [^. After his discovery, many mathematicians 
have studied this function: Hardy , @| studied this function from his 
viewpoint of the theory of elliptic functions, and Shintani pO[,pT| 



ap- 
plied it to the study on the Kronecker limit formula for zeta functions 
attached to certain algebraic fields. 

In the end of 70's, Vigneras ||2^ redefined the multiple gamma func- 
tion to be a function satisfying the generalized Bohr-Morellup theorem. 
Furthermore, Vigneras 1^^, Voros [^], Vardi and Kurokawa |I2| , 
131, |]T^, [|1^] showed that it plays an essential role to express gamma 



factors of the Selberg zeta functions of compact Riemann surfaces and 
the determinants of the Laplacians on some Riemannian manifolds. 

As we can see from these studies, the multiple gamma functions are 
fundamental for the analytic number theory: See also lOT, |]T^. How- 
ever we do not think that the theory of the multiple gamma functions 
has been fully explored. 

On the other hand, the second author of this paper introduced a 
g-analogue of the Vigneras multiple gamma functions and showed it 
to be characterized by a g-analogue of the generalized Bohr-Morellup 
theorem 
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In this paper, we will establish an asymptotic expansion formula {the 
higher Stirling formula) and an infinite product representation {the 
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Weierstrass product representation) of the Vigneras multiple gamma 
functions by considering the classical limit of the multiple g-gamma 
functions. In order to get these results, we will use the method devel- 
oped in Namely, by making use of the Euler-MacLaurin summa- 
tion formula, we derive the Euler-MacLaurin expansion of the multiple 
g-gamma functions. Taking the classical limit, we are led to the Euler- 
MacLaurin expansion of the Vigneras multiple gamma functions. The 
higher Stirling formula and the Weierstrass product representation fol- 
low from this expansion formula 

This paper is organized as follows. In Section 2, we give a survey 
of the multiple gamma functions and its g-analogue. In Section 3, we 
derive the Euler-MacLaurin expansion of the multiple g-gamma func- 
tions by using the Euler-MacLaurin summation formula. In Section 
4, we consider the classical limit of the multiple g-gamma functions 
rigorously and give an asymptotic expansion formula of the Vigneras 
multiple gamma functions. In Section 5, the Weierstrass product rep- 
resentation of this function is derived. 

Acknowledgement. The first author is partially supported by Grant- 
in-Aid for Scientific Research on Priority Area 231 "Infinite Analysis" 
and by Waseda University Grant for Special Research Project 95A-257. 

2. A SURVEY OF THE MULTIPLE GAMMA FUNCTION AND THE 
MULTIPLE g-GAMMA FUNCTION 

2.1. The gamma function. The following are well-known facts in the 
classical analysis: The Bohr-Morellup theorem says that the gamma 
function T{z) is characterized by the three conditions, 

(1) Tiz + l) = zT{z), 

(2) r(i) = i, 

(3) — logr(z+l)>0 for z>0. 
dz^ 

The gamma function is meromorphic on C, and has an infinite product 
representation 

r(.+ l)=e--^n + eH, (2.1) 

n=l ^ 

where 7 is the Euler constant. This is usually called the Weierstrass 
product formula. 
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Another representation of the gamma function is derived from the 
Hurwitz zeta function: 

^^>1- (2-2) 



k=0 

It is well-known that 

T(z) 



27r 



where ('{0, z) = f^ds, z)\s=o. 

The gamma function has an asymptotic expansion formula, i. e. the 

Stirling formula, 

logTiz + 1) ^ + log(z + 1) - + 1) - C'(0) 



+ ^[2r]2(z + ir-i' 



as 2; — >• oo in a sector A5 := G C|| argz| < n — 6} {0 < 6 < tt), 
where 

- 1 ^ n! 

n=0 

= Bk{0) (the Bernoulli number), ({s) is the Riemann zeta func- 
tion, = ^Cl-s) and [xj^ = x(x — 1) ■ ■ ■ (x — r + 1). Note that 
C'(0) = -logv4F. 

2.2. The Barnes G-function. Barnes introduced the function 
G{z) which satisfies 

(1) Giz + l) = T{z)Giz), 

(2) G(l) = l, 

(3) — logG(z+l)>0 for z>0, 
az'^ 

and he called this the "G-function" . He proved that the G-function 
has an infinite product representation. 

+ 1) = ,-.c-(o.-4.-^ n { (1 + -p + 1^) } 

and an asymptotic expansion 

logG(^ + 1) ^ (^^ _ 1^ log(^ + 1) _ _ I ^ 1 ^ ^^,(0) - log A + Oi- 

(2.0 



as 2; oo in the sector A^, where A is called the Kinkelin constant. 
Voros showed this constant can be written with the first derivative of 
the Riemann zeta function (cf ||2^, p3| ) 

logA = -C'(-l) + ^. 

2.3. The Barnes multiple gamma function. We assume that uJi,U2, ■ ■ ■ ,uJr, 
lie on the same side of some straight line through the origin on the com- 
plex plane. The Barnes zeta function |^ is defined as 

Cn{s,z;uj) ■= V — — ■ — -. 

k,M^k„=o + ^1^1 + ■ ■ ■ + 

where lj := {uji, iU2,--- , uJn). 

This is a generalization of the Hurwitz zeta function. As a general- 
ization of the formula (|2.2|) , Barnes p introduced his multiple gamma 
functions through 



exp(C(0,2;;a;)). 



where 

logp„(a;) := - lim [^(0, 2;; a;) + logz] . 

2— +0 

It is easy to see that Tn{z, uj) satisfies the functional relation 

Tn{z,Li!) _ Pn.^i(u?(z)) 

where uj{i) := {ui, ■ ■ ■ cUj-i, cUi+i, ■ ■ ■ ,uJn). 

2.4. The Vigneras multiple gamma function. As a generalization 
of the gamma function and the G-function, Vigneras |2^ introduced a 
hierarchy of functions which she called "the multiple gamma functions" . 

Theorem 2.1. There exists a unique hierarchy of functions which sat- 
isfy 

(1) Gn{z + l) = G^.^iz)Gn{z), 

(2) G„(l) = l, 

(3) ^^logG„(2;+l) >0 for z>0, 

(4) Go{z) = z. 
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Applying Dufresnoy and Pisot's results 0, she showed that these 
functions satisfying the above properties are uniquely determined and 
that Gn{z + 1) has an infinite product representation 

n-l , N 

-zE^{l) + E ^ (^i'-i(O) - E'n\l)] (2.4) 



GJz+l] 



exp 



X 

where 

EJz) 



n 



meN"-i xN* 



s(m) 



exp 



n-l 



-1) 



n—l 



n — I 



s m 



11— I 



n-l 

E 

1=0 

i)n-l{z) ■■= \0gGn-l{z + 1), 



meN"-ixN* 



-1 



\n—l 



n — I 



s m 



n—l ' 



+ (-iriog 1 + 



s m 



fc=0 



s(m) := mi + m2 + ■ ■ ■ + m„ for m = (mi, m2, ■ 

and N* = N - {0}. 

The Vigneras multiple gamma function can be regarded as a special 
case of the Barnes multiple gamma function. Namely 

Gn{z) = r„(z, (1, 1 ■ ■ ■ , I))*- ""^^ X (the normalization factor). 

In this paper, we will use the word " the multiple gamma function" to 
refer the Vigneras multiple gamma function. 

2.5. The g-gamma function. Throughout this paper, we suppose 
< g < 1. A g-analogue of the gamma function was introduced by 
Jackson [Hi , iTDI . 



r(z + l;g) = (l-g)-^n(- 



1 - 



k\ -1 



Askey |^ pointed out that this function satisfies a g-analogue of the 
Bohr-Morellup theorem. Namely, T{z\q) satisfies 

(1) T{z + l-q) = [z\,T{z-q), 

(2) r(l;g) = l, 

(3) — logr(2 + l;g) > for ^ > 0, 
where [z\q := (1 — q^)/ (1 — q). 



As q tends to 1 — 0, T{z\ q) converges T{z) uniformly with respect to 



z. A rigorous proof of this fact was given by Koornwinder |11 



Inspired by Moak's works [|T8[, the authors ||22[ derived a represen- 
tation of the g-gamma function 



1 /l-g"' 

^ogT{z] q) = {z- -) log ( ^ _ ^ 



+ log? / ^ 



^ m 

+ CM + ^^ogq + Y, 

k=l 



B. 



2k 



i2k)\ \q 



\ogq 



l-q^ 

2k-l 



(2.5) 



M2fc_i(g^ 



+ R2m{z; q), 



where 



tviU/ = logo h 

12 12g-l 



52 (t) / logg 



q'+^dt, 



R2m{z] q) 



B2rn{t) ( logg 



2m 



M2m{q'^'')dz 



(2m)! Vg*+^-l^ 

Bnif) '■= Bnit — [t]) ( [t] denotes the Gauss symbol i.e. the integral part of t), 



and the polynomial M„(x) is defined by the recurrence relation 

Mi(x) = 1, (x^ - x)-^Mnix) + nxMn{x) = M„+i(x). 

ax 



(2.6) 



Each term of the formula (2.5) converges uniformly as g ^ 1 — 0. So we 
get another proof of the uniformity of the classical limit of logr(2;; g). 

2.6. The multiple g-gamma function. Recently, one of the authors 
[|T9| constructed the function Gn{z;q) which satisfies a g-analogue of 
the generalized Bohr-Morellup theorem: 

Theorem 2.2. There exists a unique hierarchy of functions which sat- 
isfy 

(1) Gn{z + l;q) = Gn-iiz;q)Gniz;q), 

(2) G„(l;g) = l, 



(3) 



(4) Goiz;q) = [z], 



logG„+i(2 + l;g) > for z>0, 
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We call it "the multiple g-gamma function". It is given by the fol- 
lowing infinite product representation ||19| 

G„iz + 1; q) := (1 - qyU J] " ^')'"^'''^ \ 

fc=i I ^ ^ ^ / j|2.7) 

for n > 1, where 

^"("'^)=(n-l)-(n~-l)- 

In the next section, we will derive a representation of the multiple 
g-gamma function like (2.5) and consider its classical limit. This limit 
formula gives some important properties of the multiple gamma func- 
tions. 



3. The Euler-MacLaurin expansion of Gn{z + l]q) 
By means of the Euler-MacLaurin summation formula 

E /(r) = / fm + J] If {f^'-'\N) - f^'-'\M)} 

r=M -^^^ k=l 

+ (for /eC"[M,iV]), 

J M 

we give an expansion formula of the multiple g-gamma functions 
which we call the Euler-MacLaurin expansion ||2^ . This formula plays 
an important role in the following sections. 

Proposition 3.1. Suppose '^z > —1 and m > n, then 

logG„(z + l;g) 




71—1 m 

+ Y '^n,j{z)Cj{q) + Y -^Fn,r-i{z; q) - Rn,m{z] q) , 

j=0 r=l 
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where 

Fn,r-i{z;q) :-- 



n+l 



^ Y'^n /-oo 



+ 



/j+i,r--i(?) 



(n + l)! A 



1-g 



(it, 



dr-i 1 V 1 - o 



t=i 



■1) 



m— 1 /"oo 



Bm{t) 



dt"" Vn - 1 



log 



1-g 



2+1 



T/ie polynomial Gnj{z) is introduced through 



z — u 
n 



n—l 



3=0 



Proof. From (2.7) and the definition of Gn,j{z), we obtain 

logG„(;. + l-q) = -(''\ log(l - 9) - E f ~M log(l - q'^') 



(3.1) 



n— 1 



+ J]G„,(z){^fc^log(l-g^')}. 

j=0 



fe=i 



Applying tlie Euler-MacLaurin summation formula, we have 



k=l ^ ^ 



(3.2) 



n — l 



log(l - ?^+*)dt 



f /ci\'-' / -t 



n-1 



TT 1 \ Jt 



r=l 



log(l - g^+1) 



t=l 



- -^Fn,r-i{z; q) + Rn,m{z] q). 



r=l 



Let Lr{z) and be 
Li (n^) 

Lr{z) := ■ , Li{z) := - log(l - q^), U := Lr{q), 

log q 

where Lir{z) is Euler's polylogarithm 

Lir{z) J2 

Prom 



oo u 



k=l 



log(l-g^+*) = -5] 



°° qit+z)l 



I ' 



it follows that 

COD 



~t 

n — 1 



log(l - q^+')dt 



(3.3) 



n-l 



- (n-l)! ^ JF^.^'^'^' + 
where „5j is the Stirling number of the first kind defined by 

n 

where = u{u — 1) • • • (m — n + 1). Substituting (3.3) to (3.2), we 
have 



k=l ^ ^ 



^(-l>'„_i5,^(-l)'-j! 



(n — 1)! ^ a — r) 



(3.4) 



BA ( d 



r\ Wdt 



r-\ 



-t 

n — 1 



t=i 



B. 



- Y -^Fn,r-\{z; q) + Rn,m{z] q). 



r=l 



Similarly, 

oo 

^Flog(l-g^) (3.5) 

k=l 

Lr{z) and L,. cause divergence as g ^ 1 — 0, but we can prove that 
these divergent terms vanish. In order to simplify such terms, let us 
express Lr{z) as the sum of Lj. and convergent terms. 

Lemma 3.2. 

f \ — 

Li+i{z) = — lor 



l\ ^\l-q 



l-r ' ( T\lJ-r 



where 

ii r! 1 - g« 
Proof. By partial integration, we have 

Ji 1-?^ (n-r)!log g 

Let be 

- Lik+i{q) 

^k+i{z) — ^ . 

log q 

Then ( p.6|) reads 

= - E ^ (3.7) 

"^-^ / -Wn—kt'n—k iN 



fc=0 



In order to write Cn+i{z) with L^+i and T^^z), we use induction on 
n. It can be seen that 



^"^^^^^ - + ^ (n-A:)! ^^"^^ + ^ (n - A:)! ^'^^^^ 



k=0 
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which is equivalent to the claim of Lemma |3T 



Applying (3.4), (3.5) and Lemma 3.2 to (3.1), we obtain 



\ogGn{z + l;q) 



_ 7^1 



J-r 



n-1 



tr(^--o't^(^-^)' 



j=0 



r=0 



(-l)^j! 
(j-r)! 



+2 



+ 



Br f d 



r=l 



n — 1 



t=i 



n-l 



-E 



(-iy,.-iSjV-(-l)'i! 



n-l 



U " ^ ^! (j + l-r)! I 



n-l 



1-1)^-1^, ^ (-l)'j! (^+1)'+^ 
Z^ (n-l)! ^(j-/)! (/ + 1) 



n-l 

+E 

j=0 



Bj+i ( d 



{j + 1)! \dt J \n-l 



X loe 



t=i . 



1-q 



2 + 1 



1-q 



^ (-ly n-l^, (-l)^j! (-l)-(^ + 1)''V . ^ 

^ (--1)! t^O--0!t^ (l-rV. ''^'^'^ 



n-l 



r=0 ■ j=l 
- Rn,m{z; q). 
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We prove that the coefficients of the divergent terms in (3.8) vanish. 
By the following lemma, wc can see that the first term and the second 
term in (3.8) are canceled out. 



Lemma 3.3. 



^ (n-1)! Z^(j-0!^(/-r)! 

j=0 r=0 



Proof. We consider the exponential generating functions of the coeffi- 
cients of in both sides. The generating function of the left hand 
side is 

-l)^-„_i5, ^(-l)'j! ' 



(n-1)! ^ (j -/)! ^ (/-r)! r! 



z ^ u — \ 
n — 1 



On the other hand, the generating function of the right hand side is 

j=0 r=0 '' ' 

n-1 

_ (z^u- \\ 



Therefore, the coefficients of in both sides coincide. Thus, the 
claim follows. □ 
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Next, we prove the coefficient of Li vanishes. Since 

't-1 



n-l 



dt 



and 



n-l j J-. ., 

j=Q r=0 



r)\ 



n—l J-, n—1 

JD, 



r=0 j=r 



(j + l-r)! 



r=0 j=r I ^ ^ ) 



t=l 



D f / J \ ''-1 / J- 

EBr \ I a \ I z — t 



r=Q 

n-l 



dt 



n — l 



r=Q 



d\^ ^ ( z — \ 



dz 



n—l 



we have 



(the coefficient of Li) 



n 



t- 1 



n 



Br f d \ ^ ^ / z — 1 



r=0 



dz 



n — l 



-E 



n — l 



r=0 



r\ \ dz 



z=0 



We prove that the right hand side of the above formula is equal to 
zero. In a formal sense, 



__d_ \-i f d\ ^ v-^ -Br / d 



dz 



r=l 



r\ V dz 



r-l 
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Imposing the boundary condition at 2; = 0, we make the both sides 
above act on Then, 



e — 1 



-Wz-l 
— 1 



t-1 



n — 



Br 



r=l 



d 



r-1 



t - 1 

n — 1 



because (*_\) is a polynomial of (n — l)-degree. Since F{z) := — (^) 
satisfies 



F(0) = 
it can be seen that 



z-1 
n — 1 



F{z-1)-F{z) 



Thus we have the formula 



t - 1 



n — 1 



Br 



r=l 



d 

dt 



r-1 



t - 1 

n — 1 



-I t = z 



t=0 



(3.9) 



which shows that the coefficient of Li is equal to zero. Hence we have 
proved that the all coefficients of Lr vanish in (3.8). 

Finally, we calculate the coefficients of log (^ ^~i_g j ^-iid Tr{z + 1). 

Using the formula (|3l9|), we have 



the coefficient of log 



z+1 



n 



l-q 



in (3.1) 

Br 



(3.10) 



r=l 



■r-1 



Z 

n — 1 



In order to calculate the coefficients of Tr{z + 1) in (3.8), we note 
that 



z ~ u 
n — 1 



n-l 



j=0 
n-l 



(n-l)! 



^ (-1)^- r.-iS, ^ (-l)'j! >A (-l)-(z + 1)'- 
^ (n-l)! ^(j-/)!^ {l-r)\ rV 
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t=z 



i=0 



Using this identity, we have 

(the coefficient of Tr{z + 1) in (3.8)) (3.11) 



a \ I z — u 



dt J — 1 



^ \ '•-1 / 
a \ I z 



dz J \n — 1 



u=0 

Substituting (3.10) and (3.11) in (3.8), we obtain Proposition 1.1. □ 
4. The classical limit of Gn{z + 1; q) and the asymptotic 

EXPANSION OF Gn{z + 1) 

In this section, we study the classical limit of Gn{z + 1; q) using the 
Euler-MacLaurin expansion. We will see that this hmit formula gives 
an asymptotic expansion for the multiple gamma functions, which is a 
generalization of the Stirling formula for the gamma function. 

4.1. The classical limit of G'„(z+1; g). First we consiser the classical 
limit of the Euler-MacLaurin expansion in the domain {z e C\^z > 

-1}- 

Proposition 4.1. Suppose ^z > —1 and m > n. 

lim log Gniz + l;q) 

q— >1— 

C'HM-if'U}''''-" 

n-l 



R 

j=0 r=l 



where 

n+l T-, / 1 \ r-1 



r=l ^ / 



F 1 •= 

n,r— 1 • 



r-1 



n — IJ \z + 1 



t=i 



n — I 



log 



dt. 



Furthermore this convergence is uniform on any compact set in {z G 
C\^z > -1}. 



Proof. Taking Proposition |3.1j into acount, we must show that 

lim log f ^ ^ = log + 1) , (4.1) 



9^1-0^1 rl J I r! r!r ^^2) 

lim Fn,r-i{z;q) = Fn,r-i{z), (4.3) 

(jr^l — 

lim Rn.m{z]q) = Rn,m{z), (4.4) 

lim Q(g) = Q, (4.5) 

q— >1— 

and further have to show that this convergence is uniform. Here we 
prove only (4.4). The other formulas can be verified in a similar way. 
Since 

.^^"-o (I) { G - 1) (r^) } = (I) { (n - 1) (fri 

in order to show (4.4), it is sufficient to prove that the procedure of tak- 
ing the classical limit commutes with the integration. Let us introduce 
polynomials My.{x) through 

^ log(l - g^+*) = - (y^) ' q'^'MM'"-') 

(cf.[18]). They satisfy the recurrence relation. 

Mi{x) = 1, (x^ - x)-^M„(x) + {(r - l)x + l}M^(x) = M^+i(x) 

ax 

and Mr{l) = (r — 1)!. Using this we have 
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(n 



^|{m)'..-.4(T)w. 



X / BUt)t^~' (r^) I'^'Mm-iiq'^ndt 



1 

Therefore we have to show 

m—l 



Lemma 4.2. Suppose that a G N 5e fixed and that i/q and yi are fixed 
constants such that —1 < uq < yi. Then there exists a constant C 
depending on and yi such that 

log g \ " C 



< 



qt+y _ 1 



< 



{t + yY 



for y G [yo, 1 < t < oo, < g < 1. 



Proof. Letting q = e we define 



{t+y)5 



-s{t+y)-i 
1 



-5{t+y) 



ipit,y,6) := 

Then, ip(t,y,S) is positive, continuous and bounded on 



6>0 
6 = 0. 



D :-- 



it,y,6) 



l<t< +00, 6>0, yo<y<yi}. 



So we can put C := supip{t,y,6). 

D 



□ 



Suppose — 1 < 2/0 < < yi. By Lemma [4 .21 , there exists a constant 
C depending only on yQ and yi such that 

m— i 



1 - 



rn—l 



for 1 < t < +00, < g < 1. 

Since m — j > 2, the right hand side above is integrable over 1 < 
t < +00. Therefore, Lebesgue's convergence theorem ensures (4.6). 
We have thus proved the limit formula (4.4). Next we show that this 
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convergence is uniform on any compact set in the domain {z & C\^z > 

-!}• 
Put 

Prom the consideration above, we see that there exists a constant C 
depending only on yo and yi such that 

/■°°_ dt 

mz,q)\<Cj B^{t)- for ?lizE[yo,yi] and < g < 1. 

Hence {^{z,q)\0 < q < 1} is a uniformly bounded family of functions 
and 



/CO 
Bm{t) 



— , A^-\-\r-\m-l-\)\ 



dt as g — > 1 — 0. 



By Vitali's convergence theorem, this convergence is uniform on any 
compact set in the domain {z G C|3fJ2; > — 1}. □ 

The constant Cj in Proposition 4.1 can be expressed in terms of the 
special value of the Riemann zeta function. 

Lemma 4.3. 

1 



Cj = -C'i-j) 



Proof. Prom the definition of Cl^), 

C'W = -f;i^ tor &>1. 

fe=l 

By the Euler-MacLaurin summation formula, we obtain 

1 r°°— f dV 



(4.7) 



Since 



7 \ »l r\ 
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(4.7) can be analytically continued to {z e C\'^z > —n + 1}. So if we 
put s — —j, n — j + 2, then we obtain 



r=l 



t=l 



(i + i)^ 



□ 



Next we prove that the limit function in Proposition 4.1 coincides 
with the multiple gamma funciton. 

Theorem 4.4. Suppose m > n. Then, as q ^ 1 — 0, Gn{z + 
converges to Gn{z + 1) uniformly on any compact set in the domain 
C\Z<o and 



logG„(^ + l) 



(4.8) 



z + 1 
n 



r=l 

n-1 



E 

r=l 

r-1 



Br. f d 



r! V dz 



r-l 



log(. + l) 



dz J \n — lj\ r\r 



- EGn,(.) |C'(-J) + I + E ^i^n.-l(.) 



Rn,m (^) • 

Proof. In the domain {2; G C|3?2; > —1}, we have already proved the 
existence of the limit function and uniformity of the convergence. Let 
us put 

Gn{z + l) := lim G,,{z + l]q). 
g— >i— 

Because of the uniformity of the convergence, we have particularly 



lim 

g-»l-0 



n+1 



{logG„(;^ + l;g)} 
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7 \ "+1 



so that, from Theorem 2.2, Gn{z+1) satisfies the conditions in Theorem 
2.1. Namely 

(1) Gn{z + 1) = Gn^l{z)Gn{z) 

(2) — logGniz + 1) > for z>0 



(3) G„(l) = l 

(4) G^{z + l)=z + l. 

Since a hierarchy of such functions is uniquely determined, so Gn+i{z + 
1) = Gn{z + 1) in {z G C|3fJ2; > —1}. Thus the claim of the theorem 
in the case that "Rz > — 1 has been proved for {z E C\^z > — 1}. 
Next, we show that in {z G C|3?2; < —1, z ^ —1}, 

Gn{z + l]q) ^Gn{z + l) as 

and that the convergence is uniform on any compact set in this domain. 
For the proof, we use induction on n. 

The case that n = 1 was considered by Koornwinder |[Tl|. Let K be 
a compact set in {2; G C| — 2 < < —1, z 7^ —1}- If q is sufficiently 
close to 1 then [z + 1]^ 7^ on K. Therefore as g ^ 1 — 0, 

T(z + 2-,q) 

uniformly converges to 

T{z + 2) 
z + l 

on K. 

We assume that 



r(^ + i) 



Gn-i{z + 1; g) — s> Gn-i{z + 1) as g — 1 — 

and that the convergence is uniform on K 

From (2.7), we see that, if q is suffciently close to 1, Gn^i{z + 1; q) 
has no poles and no zeros on K, neither has Gn-i{z + 1) from (2.5). 
Therefore, as g — > 1 — 0, 

^"^" + '''^-G._,(. + l;g) 
uniformly converges to 

G„_,(. + l)-^"(" + ') 

on K. 
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Repeating this procedure, we can verify, for any n, Gn{z + l;q) con- 
verge to Gn{z + 1) in a compact set in the domain {—3 < < —2, z ^ 
—2}, {—4 < < —3, z 7^ —3}, • • • . Thus the claim of the theorem is 
proved. □ 

4.2. Asymptotic expansion of Gn{z + 1). Let us call the expression 
(4.8) the Euler-MacLaurin expansion of Gn{z + 1). We should note 
that (4.8) is valid for z G C\{^z < —1}. We show that it gives an 
asymptotic expansion of Gn{z + 1) as \z\ — >• oo, i.e. the higher Stirling 
formula. 

Theorem 4.5. Let < 6 < tt, then 

l0gGn{z+l) 

■r).|:5B)"(.i.)}««. 

n-1 



r 1 1 R 

as \z\ oo in the sector {z e C||arg'2;| < tt — S}. 
Proof. Straightforward calculation shows that 



|)'log(. + *) 

t=l V. " ^ J t=i 



1=1 \k=0 



{z + iy 



Thus, if r > n, then Fn^r-i{z) = 0{z ^~^^) as \z\ — >• oo. So we can 
see that 



-^n,2r-3(^) 

Furthermore we can see that 
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-1 /■°°_ fd 



2m 



2m-j 



Noting that 



1 1 

< 



\z + t\ |t||sin(5| 



in the sector {z G C||arg'2;| < tt — (5} and that \B2m{t)\ < |-B2m| for 
< t < 1, we have 



2m 



(2m)! 



I J |sin(5|2™-J i2m-j 



Hence, 
as l^l ^ oo in the sector. 



□ 



Let us exhibit some examples of the higher Stirhng formula. In the 
case that n = 1, we obtain 

logG'i(;2 + l) =logr(;2 + l) 

~ (^z + log(^ + 1) _ + 1) _ c'(0) 



B 



2r 



- [2r]2 {z + 1)2-1 



This is the Stirling formula since ('{0) — log(27r). 
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Furthermore, in the case that n = 2, we obtain 



logG2{z + l) 

log(2; + 1) h - 

2 12/ 'a 2 4 

-<'(o) + c'(-i) 

" 12 ITT + ^ [2r]3 + 1)2-1 - 2r + 1), 



which coincides with the formula (2.3). In the case that n = 3,4, and 
5, we have the following results. 

Proposition 4.6. The higher Stiring formula for n = 3, 4 and 5 are 
as follows: 

\ogGs{z + l) 

'z^ z^ 1 \ , / .X 11 s So -2 13 
6-4+24j^°S^^ + ^)-36^ +24^ +3-72 

- ^C'(o) + ^C'(-i) - ^C'(-2) 

+ + Y {z^ - (6r - ll)z + (4r2 - 16r + 16)| 

12^ + 1 ^\MAz^\Y'-^^ \ ) \ )f 



r=2 



\ogG^{z + l) 

' 19 



24 6 6 720 



log(^ + 1) 



4 4 2 3 ^2 11 31 

— Z^ + 77^ + ^ - ^Z + 



72 9 8 36 144 

£!^i£!±^C'(o) + 5£!z|i±^c'(-i) - ^C'(-2) + ic'(-3) 

11 1 1 / 2 13 5 



12^+ 1 720(^ + 1)3 V '2 2 
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°° B 1 

+ E [2^ (^ + 1)2.-1 - (12r - 27)^2 + (20r2 - 94r + lll)z 
-(»r" - 56r2 + 134r - 109)} . 



r=3 

,3 



logG'5(^ + l) 



137 , 39 4 461 o ^2 323 5639 



7200 320 2160 1440 1440 43200 
- 6^=^ + 11^2 _ 4^3-18^2 + 22^-6 , 

24 ^ ^ 24 ^ ^~ 

- ""-^f + ^ C-M) + ^C'(-3) - ic'(-4) 

1 1 1 1 /35 2 45 9\ 

^ 12^ + 1 " 720 (^ + 1)3 VT^ ^T^^2y 

oo ^ 1 

+ V 1-2^ - (20r - 54) ^3 + (70r2 - 375r + 506) 



200 o , 4420 

— - 540r2 + — — r - 1354 z 



4 536 o 9 4279 
+16r^ - — + 754r2 —r + 1021 



5. The Weierstrass product repesentation for the 

Gn{z+1) 

By calculating the formula (4.8) in the case of m = n + 1, we de- 
rive the Weierstrass product representation for the multiple gamma 
functions. Main theorem of this section is the following: 

Theorem 5.1. For n e N, we have 

+ 1) = exp {F^{z)) n Ul + l) exp k)) , 

k=l ^ ^ 
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where 



n—1 



n-2 



3=0 



r=0 



1 f^-'' 

r\ \du j \n — 1 



X C'(-r) 



u=0 



z — u 
n — 1 



du X 7, 



^ n— 2 ^' n—1 



i+1 



r-=0 



r=l ^ 



i+1 



B, 



r=0 



-r+1 



logt 



dtj [j + l ° + 



t=i 



Proof of this theorem will be carried out through the sections 5.1 ~ 
5.4, and some examples of this representation will be discussed in the 
section 5.5. 

5.1. Rewriting the Euler-MacLaurin expansion of G„ (2; + 1). In 

this section we prove the following proposition. 



Proposition 5.2. 



n-l 



j=0 

where 

Kjiz) log(. + 1) - C'(-J) + Pjiz + 1) 



(5.1) 



J + 



k=l 



Pj{z + k + l) - Pj{z + k) + 



B,+i(-- + /,- + l), fz + k + 1 



log 



z + k 
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Furthermore the infinite sum of the last term is absolutely convergent. 



Proof. Prom the definiftion of Gn^ri^), we have 

r— 1 / \ / 7 \ r—l 



d 

di 

and hence 



z 

n — 1 



E - 



r=l 



d_ 

dz 



r-l 



d_ 

du 



z 

n — 1 



z — u 
n — 1 



{r-iy.Gn,r-liz) 



u=0 



x—{iz + iy-i} 

r ir 



n-l 



j=0 



z=0 

{Z + 1)-^'+^ - 1 



(5.2) 



(5.3) 



Since 



n-l 



z-{z + t) 
n-l 



n—l 
j=0 



we have, for 1 < r < n — 1 



n— 1 



j=0 



(5.4) 



In order to calculate the coefficient of log(2; + 1), we use the following 
lemma. 



Lemma 5.3. 

fz + l 
\ n 



n— 1 



n-l 



3=0 



j=0 



Proof. Since the both sides are polynomials in z, it is sufficient to prove 
that the formula holds for a sufficiently large, arbitrary integer N. 
Noting that 

Bj^,{N + 1) - Bj^, _ [Ez=i^^' (J>0) 
j + 1 [N+l (j = 0), 

we have 

EG„.w5^±i^^-EG„,(7V)^ 

26 



N n-1 

5^5^G„,(iV)/^+G„,o(iV) 

1=1 j=0 



N 

E 

1=1 



72—1/ \n — 1 



+ 



n I \n — 1 

N + 1 
n 



□ 



From (|5.2|) and Lemma 5.3, it follows that 



n-l 



/ j=0 ■' ^ 



(5.5) 



Next we calculate Rn.n+i{z)- From the definition of Gn.r{z)., we have 

Rn,n+l{z) 



(n + 1)! A ^"+^^^HI 



n+l 



z - (2; + t)\ + t 



n-l 
3=0 



n — 1 

n+l 



log 



z + 1 



(2 + ty log 



n-l 



j=0 



n+l 



+ 1 + ky log 



^ + 1 . 
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Here we have, for j < n — 1, 

n+l 



(5.7) 



so that the infinite sum in (5.6) converges absolutely because j—n—1 < 
-2. 

By means of the Euler-MacLaurin summation formula, we have 



(. I! (I) { ^ ^ ^ (^) } 



(5.8) 



^ ^ + E ^b1.-i(^ + k + 1)^+^-^ [ log(;. + k+l) 



+ 



(z + ky + 



r=l 

{z + k + iy+^ _ {z + ky+^ 

JTl 3+1 



1 



log(2; + 1) 



- E ^'^^•'^ {(^ + ^ + - (-^ + A;)^'+'-''} . 



r=0 ■ 



The coefficient of log(2; + 1) of the above formula vanishes because 
of the following lemma. 

Lemma 5.4. 

7 + 1 7 + 1 ^ ' 



i+1 



r=l 
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Proof. The first equality follows from the identity 



r=0 



J + 1 



+ k + ly+'-'Br. 



The second equality holds because of the Euler-MacLaurin summation 
formula for {z + t + ky on [0, 1]. □ 



Applying Lemma 5.4 to (5.8), we obtain 

71—1 OD 

Rn,n+l{z) = — Gn.j{z) 
j=0 



k=l 



Bj+i{z + k + l)^ fz + k + 1 
— log ' 



J + 1 



z + k 



n+l 



B. 



+ E TT ^^vr {{z + k + ly^'-^ -{z + ky^'-^} 



r=0 



If r > j + 2, then {z + k + -{z + ky^^'"' decreases more rapidly 

than k^^ as A; — > cxo. Thus we have 



n— 1 
j=0 



n+l „ 



r=j+2 



Tl 



(5.9) 



k=l 



Bi+i{z + k + l)^ fz + k+l 
- log ' 



J + 1 



z + k 



Pj{z + k + l) - Pj{z + k) 



and the infinite sum is absolutely convergent. 
Substituting (5.3), (U), (5.5), (5.9) to (4.8), and noting that 



, n+l „ n+l „ 



(j + 1)^ ^ rl 
we obtain (5.1). 



r=j+2 



□ 



The next lemma which will be useful in the section 5.3 is deduced 
from the above proof. 

Lemma 5.5. 

5,+i(z + fc + 1) , fz + k + 1 
- log ' 



J + 1 



(^^±^^+P,iz + k + l)-P,iz + k) 



= 0{k ) as k oo. 

In other words, the polynomial Pj{z + k + 1) — Pj{z + k) is a convergent 
factor such that the infinite sum in (5.9) converges absolutely. 
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5.2. An infinite product representation for the C{~j)- We de- 
rive an infinite product representation for C{~j) the same way as 
in the section 3.1. 



Proposition 5.6. 

exp(C'(-j)) = exp(P,(l)) n U 1 + ^] exp {P,{k + 1) - P,{k)) 



and the infinite product converges absolutely. 



Proof. From the proof of Lemma 4.3, we have 



U + 2)! A 



(j + 2)! A '^'''\dtj 
Furthermore we can see by the same way as (5.7) that 

J Bj+2{t)i-^] {tnogt}dt = 0{k-^) as k-^oo, 
Hence we have 



(5.10) 



(j + 2)!7i 



BH2it)[j^j {tnogtjdt 
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and the infinite sum of tfiis formula is absolutely convergent. In the 
same fashion as in the previous section, we have 



J + 



(j + 2)! V't ^+1 
From (|5.10 ) and (5.11), it follows that 
S,+i(fc + l) 



3 + 



log (l + ^) + ^.(^ + 1) - ^.(^) (5.12) 



0{k ^) as A; ^ oo. 



Hence we have 



c'(-j) = p,(i) + Ei°g 
fe=i 



1 + ^) exp{P,(fc + l)-P,(fc)} 



and the infinite sum is absolutely convergent. 



□ 



Remark 5.7. In a similar fashion to Proposition 5.6, an infinite prod- 
uct representaion of C'{—j,z) (:= -^((s, z)\s=-j) can be given as fol- 



lows: 



exp(C'(-j, z)) = exp 



J + 1 



\ogz + P(z] 



X 



n 

fc=0 



-Bj + l(z+fe+l) 

Z + + 1\ 



z + A; 



exp {Pj{z + k + 1) - Pj{z + k)) 



5.3. Good representation for Kj{z). In this section, we give a "good" 
representation for Kj{z). The following two lemmas are useful in sub- 
sequent arguments. 
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Lemma 5.8. Let k be a positive integer and k —>■ oo, then we have 

'^y^ ur-l _ miu-l 

1=1 

j+i r+l 



(1) k^\og{l + ^)+J2^-^k-^ = Oik 

1=1 

Bi.Uz + k + l)^ / z\ 1 4^ .sv^(-l)'. 



0{k 



r=0 1=1 

-2\ 



Proof. (1) is clear. (2) follows from the identity 

i+i 

B,+,{z + k + l) = J2 B,+i-r{l){z + ky. 

T=0 



Lemma 5.9. There exists a unique polynomial A{k, z) such that 

5!±l(l±i±il.og(l + i)+^(fc..) = 0(t-), 

A{k, z) is equal to Pj{z + k + 1) — Pj{z + k) . 
Proof. By integrating the both side of 



□ 



i+i 

i?,+i(^+l) = ^5,+i_.(l) 

r=0 

from —1 to 0, we get the formula 

i+i 



Prom Lemma |5.8| (2) and (|5.13|) , a polynomial A{k, z) satisfying 
^2±}!^±^±}llog(l + '-]+A{k,z)=0{k'^) as fc^oo, 

is uniquely determined and the polynomial Pj{z + k + 1) — Pj{z + k) 
satisfies (|5.14|) by Lemma |575| . Therefore, A{k, z) = Pj{z + k + 1) — 
Pj{z + k). □ 

Using these lemmas, we prove the following proposition. 
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Proposition 5.10. Let k be a positive integer and define 



' l,r-=0 ^ ^ 1=1 

:= -{I)k+i, 

1 ^^■+1 

(///), := J2 ^''^ {Prik + 1) - Prik)} - --— 
r=0 '^^^^ 

Then we have 



and 



j+1 \ z^k 



^^-(^ + ^) log )+(/), 



(5.15) 



^ (-lyv 

7 + 11^ r + 1 



i + i fn; '■ + 1 '+"^mU 4+1 



Furthermore each term in the right hand side decreases like 0{k ^) 
as k ^ oo. 



Proof. Making use of the identity, 

Bj+i{z + k)- Bj+i{z + k + l) 
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-{z + ky, 



we have 



j+1 \ z+k 

Bj+i{z + k)^ f k \ + + , / z + A; + 1 
log —T + ^ log 



j + 1 \z + kj j + 1 \ k + 1 

+ A; + 1) , A + lA , fz + k 

+ ^ log ( J ~ ^ 



By Lemma |5.8| (1), (5.12) and the identity 



Bj+i{z + k + l) _ 1 v"/j + l 



r.— n ^ ^ 



we can see that each term in the right hand side of (5.15) is 0(/c~^). 

On the other hand, it can be seen that {I)k ~ {IV) k are polynomials 
of z and that 

1 ^v^/i + l^o .^(-lr^^+' 



r + 1 I A; 



+ (a polynomial of k), 



+ (a polynomial of k), 

z^+'^ 1 

{III)k = T + polynomial of k), 

3 + Ik 



r + 1 A; + 1 



(^^)^ = I E 7^ (^,) I ^ + (polynomial of k) 
+ (a polynomial of A;). 



r=0 
^^■+1 1 



3 + Ik 
Hence, if we put 
B{k, z) := (/), + (//)fe + (///), + {IV)k 
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then, B{k, z) is a polynomial of fc, z and it satisfies 

:^±i(£±A±l)iog(^l + ^)+i?(A;,^) = 0(/t-2) as k^oo. 

By Lemma 5.5, we have 

B{k.z) = Pj{z + k+l)- Pj{z + k). (5.17) 
By (5.16) and (|5.17| ), we can deduce (5.15). □ 



Proposition 5.11. 



K,{z) = Q,{z)+Y,(^^z^-'X'{-r) 



r=l 



+ ^ -(-- + *)Mog (i + + 5: (^^) j: tiriik' 



k=l r=0 ^ ^ Z=l 

anc? t/ie infinite sum of this formula converges absolutely. 



Proof. The summation of each term in (5.15) from A; = 1 to A; = oo, 
is absolutely convergent because of Proposition 5.9, and the following 
caluculation is possible. Namely, 

f:l^.os(^)..)4 



k=l 



k=l 



J + 1 \ k+1 



■' r=0 ^ ^ l=\ 



Noting Theorem ^.6| and 

7 

we have 



f:{'o.(i+i) 



r=0 
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By (5.18) and (5.19), we obtain 

Bj+i{z + k + l) fz + k + 1 



E 

k=l 



J + 1 



J + 1 



log 



z + k 



+ Pjiz + k + 1) - Pjiz + k) 

{5.21 



\og{z + 1) 



r=l ^ ^ 1=1 



+ ^ p)^^-^ {C'(-r) - P.(l)} - 



r=0 



J + 1 



+ ^ ^ + log (l + i) + ^ (^J 5; 



r=0 



«=1 



The proof is completed by substituting (5.20) to the definition of 
Kj{z) in Proposition ^.2| . □ 



5.4. A proof of main theorem. Using the results in the sections 5.1 



~ 5.4, we prove Theorem |5.1| . By Proposition ^.11| , we have 
logG„(z + l) (5.21) 



n-l 



^.W + EC)-^-^C'(-r)-|9^7 



+E 

k=l 



+ *)Mog ( 1 + i-) + E (!) E ^-^^k^-' 



It is easily seen that 

n-l j-l 



k 



r=Q 



1=0 



I 



z^-^Ci-r) 



(5.22) 



r=0 
n-2 rn-l 

r=0 lj=0 
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n-2 

E 

r=0 



1 / (9 \ fz-u 
r! \du I V^^ — 1 



X C'(-r), 



?i=0 



(5.23) 



n-l 



n — 1 



(5.24) 



Substituting (5.22), (5.23), (5.24) to (5.21), we have 



(5.25) 



n-l 



n-2 



^G„,(z)Q,(z)+^ 

j=0 r=0 



1 / d \ ^ I z — u 



r\ \du J \n — 1 



X C'(-r) 



n=0 



2; — M 
n — 1 



du X 'y 



E 

fc=i 



) n-l j , r+l 

|°g(i + i) + EG,wwEC>'-E 
j=0 r=0 ^ ^ Z=l 



Thus, in order to prove Theorem |5.1| , it is sufficient to show 

n-l j 



j=0 



n— 2 f n— 1 



r=0 /=1 



n— I'S'j' 



(5.26) 



^ ' u=-l Lr=M+l 



Since (5.24) and 



+ *r log (i + 1) + 1 (;\- f ^>i^^r-' 

r=0 ^ ^ /=1 



0{k ^) as A; — > oo, 
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we have 

n-l 



^/ ^ j=o r=o V/ i (5.27) 

= 0(A;~^) as k ^ oo, 
while we obtain 

= Oi^k"^) as k ^ oo, 



by Lemma [57^ (1). We can deduce that (5.27) and (5.28) imply (5.26) 
by the same arguments as in Lemma 5.8. Hence the proof is completed. 

□ 

5.5. Examples of the Weierstrass product representation for 

Gn{z + 1). We give some examples of the Weierstrass product repre- 
sentation for the multiple gamma functions. 
In the case that n = 1, we have 

°° r -1 
G,iz + l) = Tiz+l)=e-^'l[Ul + ^y e-i 

k=l ^ 

This is the Weierstrass product representation for the gamma function. 
In the case that n = 2, we have 

G,iz + 1) = Giz + 1) = e-<'W-4--^ n I ^ 'k) + k) } 

Since C'(0) = — | log(27r), this is the Weierstrass product representation 
for the Barnes G-function 

In the case that n = 3, 4 and 5, we obtain the following results. 

Proposition 5.12. The Weierstrass product representations in the case 
that n = 3, A and 5 are as follows: 

G,{z+1) 

- {- 4 + ¥ + 24^ + ^ - ^C'(O) - - 4 7 
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X 



n 

k=l 



On) 



fc(fc+i) 

2 



exp 











4y 





' z\ z , 

4-2+2^ 



^4(^ + 1) 



61 4 13 o 19 2 5 
= ^"P^144" +18" +144" -24" 



z^,, z'^-2z^,, , z^-^z^ + 2z^,, , z^-Az^ + Az'^ 

+ ^-c'(-i) ^ — c'(o) — 7 



X 



n 

ik=l 



fc(fc+l)(fc+2) 



exp 



24 6 + 6 j A; 



18 4 sJ + VlS 6^ 



^5(^ + 1) 



= exp 



288 



, 7 . 173 , z^ 2827 
-2 + TTT^ - + z 



64 



864 



+ ^C'(-3) - ^^C'(-2) + 



36 17280 

2z^ - + llz., 



12 



C'(-i) 



- 6z=^ + 11^2 - 6^ , 6^5-45^^ + 110^=^-90^2 

C (0) ^ 7 



24 



720 



X 



n 

fe=i 



1 + 



fc(fc+l)(fc+2)(fc+3) 
^ \ 24 



exp 



^4 11 ^ ^ ^2X 1 

120 ~ 16 + 72" ~ "8 / ^ 



_ 11 2 
96 ~ 12 + 48^ 



(z^ z^ 11 
172 "8+24'^ 



24 



24 
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